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Abstract 

A family of reflected Brownian motions is used to construct Dyson's process of non-colliding 
Brownian motions. A number of explicit formulae are given, including one for the distribution 
of a family of coalescing Brownian motions. 

1 Introduction 

The ordered eigenvalues Y\ (t) < Yz (t) < . . . < Y n (t) of a Brownian motion in the space of n x n 
Hermitian matrices form a diffusion process which satisfies the stochastic differential equations, 



(1) Yi(t) =y t + fc(t) + J2 



ds 



where 0x, (3%, . . . , j3 n are independent real Brownian motions. This is a result that goes back to Dyson 
[7] and we will refer to Y as a Dyson non-colliding Brownian motion. A number of important papers 
in recent years have developed a link between random matrices and certain combinatorial models, 
involving random permutations, last passage percolation, random tilings, random growth models 
and queuing systems, see Baik, Deift and Johansson, P and Johansson, amongst many others. 
A recent survey is given by Konig J3| . At the heart of this connection lies the Robinson-Schensted- 
Knuth algorithm, a combinatorial procedure which has its origins in group representation theory, 
and using this the following remarkable formula, was observed by Gravner, Tracy and Widom, [TU] 
and Baryshnikov [2], representing the largest eigenvalue Y n (t) ( assuming Y(0) — 0) in terms of 
independent, real- valued, Brownian motions Bi, B 2l ■ ■ . , B n , 

n 

(2) Y n (t) d ^ sup r^W-Bi^i)}. 

o=t <t 1 <...<t„=t f 
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O'Connell and Yor, |TS], give a proof of this identity by considering reversibility properties of a 
queuing system, which in a subsequent paper, O'Connell 16 , is shown to be linked to the RSK 
algorithm also. Another proof, again involving RSK, is given by Doumerc, 

In this paper an a different proof of the identity @ is given, based around the following con- 
struction. Let (Y(t);t > 0) be a Dyson process, with components Yi,Y 2 , . . .Y n solving Q). Let 
(X(t);t> 0) be a process with (n + 1) components which are interlaced with those of Y, meaning 

(3) X 1 (t)<Y 1 {t)<X 2 (t)...<Y n {t)<X n+1 {t), foralH>0, 
and which satisfies the equations 

(4) X % (t) =x l+ll (t) + { L7 (t) - L+ (t) } . 

Here 0y(t);t > 0) is a standard Brownian motion in R n+1 , independent of the Brownian motion 
j3 which drives Y. The processes (Lf(t);t > 0) and (L~(t);t > 0) are continuous non-decreasing 
processes that increase only at times when Xi(t) — Yi(t) and Xi{t) = Yi-%(t) respectively: they are 
twice the semimartingale local times at zero of Xi—Yi and Xi — l^-i. The two exceptional cases L^it) 
and L^ +1 (t) are defined to be identically zero. Conditionally on Y the particles corresponding to X 
evolve as independent Brownian motions except when collisions occur with particles corresponding 
to Y. Think of the particles corresponding to the components of Y as being "heavy" so that in 
collisions with the "light" particles corresponding to components of X their motion is unaffected. 



On the other hand the light particles receive a singular drift from the collisions which maintains the 
interlacing. We will verify that is possible to start X and Y from the origin so that xi = yj = for all 
i and j. Then, see Proposition [S] the process X is distributed as a Dyson non-colliding process with 
(n + 1) particles. Thus if we observe only the particles corresponding to the components of X, the 
singular drifts that these particles experience from collisions with the unseen particles corresponding 
to Y are somewhat magically transmuted into an electrostatic repulsion. This is a consequence of a 
relationship between the semigroup of the extended process (X, Y) and the semigroup of X that is 
called an intertwining relation. 



2 A duality between interlaced Brownian motions 

Consider a continuous, adapted, R™ +1 x R™-valued process (X(t),Y(t);t > 0) having components 
Xi(t), X%(t), . . . X n +i(t) and Yy(t), ^(t), • • • ,Y n (t) which is defined on filtered probability space 
(fl, J 7 , {jF t } 4>0 , Q™ y ) satisfying, for all t > 0, the interlacing condition 

Xx(t) < Yt(t) < X 2 {t) <...< Y n (t) < X n+1 (t), 

and the equations 

(5) Y i (t) = y i +p i {tAr), 

(6) Xi(t) = Xi + ii(t At) + LT(t At) - Lf(t At), 



where, 

t is the stopping time given by r = inf{f > : Yi(t) — Y i+1 (t) for some i G {1, 2, . . . , n — 1}}, 
(3\, @2, ■ ■ - fin, 7ij 72, ■ ■ ■ ln+i are independent JF t -Brownian motions, 

Li(t) = L^ l+1 {t) = for all t > 0, otherwise the processes Lf and L~ are continuous, non-decreasing 
and increase only when Xi = Yi and Xi = Y^\ respectively, 

L+(t) - f l(Xi{s) = Y{s))dL+{s) Lr(t) = f l(Xi(s) = Y^dLr ( s ). 
Jo Jo 



The process just defined is called a stopped, semimartingale reflecting Brownian motion. For general 
results on such processes see, for example, Dai and Williams, 0]. In this case it is not difficult to 
give a pathwise construction starting from the Brownian motions f3i, for t 6 {l,2,...,n}, and 7^ for 
i G {1, 2, . . . , n + 1}, together with the choice of initial co-ordinates x\ < yx < x% . . . < y n < x„ + i. 
We obtain Yi immediately. Xi is constructed by alternately using the usual Skorokhod construction 
to push Xi up from Yi-i and down from Yi. For more details see Section 3 of ^H], where a similar 
construction is used. In fact by the same argument as Lemma 6 in Soucaliuc, Toth and Werner, |19| 
pathwise uniqueness holds, and hence the law of (X,Y) is uniquely determined. This uniqueness 
implies, by standard methods, that the process is Markovian, and in fact we are able to give an 
explicit formula for its transition probabilities. 

We denote by <p t the centered Gaussian density with variance t. <!> t is the corresponding distri- 
bution function 

Mv)= [ -^=cxp{-z 2 /(2t)}dz } 

J -00 V 27TI 

and 

Let W n+1 ' n = {(x,y) G R ,l+1 x R™ : x x < yi < x 2 < . . . < y„ < x n+1 }. Define q?({x,y), (x',y')) 
for (x, y), (x' , y') G W n+1 ' n and t > to be equal to determinant of the (2n + 1) x (2n + 1) matrix 

(A t {x,x') B t {x,y')\ 
\C t {y,x') D t (y,y')J 

where 



A t (x, x') is an (n + 1) x (n + 1) matrix with (i, j)th element <fit(x'j — Xi); 

B t (x,y') is an (n + 1) x n matrix with (i, j)th clement &t{y'j — %i) — l(j > 

Ct(y,x') is an n x (n + 1) matrix with (i,j)th element <j>' t {x'j — yi); 

D t (y,y') is an n x n matrix with (i, j)th element 4>t{y'j — yi)- 

Lemma 1. For any f : W n+1 ' n — > R which is bounded and continuous, and zero in a neighbourhood 
of the boundary ofW n+1 ' n , 



lip / q?(w,w')f(w')dw' = f( 



w), 



uniformly for all w = (x,y) G \y n + 1 , n _ 

Proposition 2. (</";£ > 0) are a family of transition densities for the process (X, Y) fciZ/ed at ifte 
instant t, that is to say, for t > and (a;, y), (x',y') G W n+1,n , 

q?((x,y),(x',y'))dx'dy' = Q^(X t G aV,Y t G dy';t < r). 

Proof. For any choice of z' G R, each of the functions (t, z) i— > $ t (z' — z), (t, z) i— » 0t(z' — z) and 
(i, z) i — ► <p' t ( z '~ z ) satisfies the heat equation on (0, oo) x R. Thus, by differentiating the determinant, 
we find that, 

2n+1 f) 2 n n F)n n 

( ? ) 2 E = ^H^') e (°'°°) x R2n+1 x R2n+1 - 

i=l Wi 

We need to identify certain boundary conditions. We treat w' = (x',y') G W n+1 ' n as fixed. First 
consider (x,y) G dW n+1 ' n satisfying ?/, = y i+ i for some « G {1, 2, ... n — 1}. We see that the ith 
and (i + l)th rows of both C t (y, x') and D t (y, y') are equal, and hence q?((x, y), (a/, y')) = 0. Next 
consider (x, y) G dW n+1 ' n satisfying Xi = yi for some i e {1,2, ... n}. Calculate ^-q"((x, y), (V, y')) 
by differentiating the ith rows of A t (x,x') and B t (x,y'). Notice that, under our assumption that 
Xi = yi, the ith row of -£^A t (x,x') is equal to the zth row of —C t (y,x'). Likewise the ith row of 
■£^B t (x,y') is equal to the ith row of —D t (y,x'). Thus we deduce that T^-g"((x,y), (x',y')) = 0. 
Finally consider (x,y) G dW n+1,n satisfying Xi+\ = yi for some i G {1,2, . . .n}. Similarly to the 
previous case we obtain d ® q™((x,y), (x',y')) = 0. 

Let / : W n+1 ' n -^Rbea bounded and continuous, and zero in a neighbourhood of the boundary. 
Then define a smooth function F on (0, oo) x W n+1 ' n via 

F(t,w)= [ q?(w,w')f(w')dw'. 

By virtue of the above observations regarding g™, and differentiating through the integral, we find 
that 

^E^ w ) = ^ ffl ) on(0,oo)xW"+ 1 '«, 
with the boundary conditions 

F(t, w) = whenever w — (x, y) satisfies yi = yi + \ 
dF 

— — (t, w) = whenever w = (x, y) satisfies Xi = y-i 

OXi 

dF 

— (t, w) = whenever w = (x, y) satisfies x i+ i = yi 

ox i+ i 

Fix T, e > 0. Applying Ito's formula, wc find that the process (F((T + e - t, (X t ,Y t )) ; t G [0, T]) is 
a local martingale, which is easily seen to be bounded and hence is a true martingale. Thus 

F(T + e, (x, y)) = Q» „ [F(e, (X T , Y T ))] 

= Ql y [F(e,(X T ,Y T ))l(T<r)]. 



Appealing to the previous lemma, we may let e J, and so obtain, 



F(T, (x, y)) = Ql y [f(X T , Y T ))1(T < r)] . 

Since the part of the distribution of (Xt, Yr) that charges the boundary of W n+1 ' n exactly corre- 
sponds to the event {T > r} this suffices to prove the proposition. □ 

We now consider a second reflected semimartingale Brownian motion (X,Yj having components 
Xi(t),X2(t), . . . X n+ i(t) and Yi(i), ^(i), • • • , Y n (t) which is defined on filtered probability space 
(f2, !F, {jFt} 4>0 , Q" y ) satisfying, for all t > 0, the interlacing condition 

Xx(t) < Yx{t) < X 2 {t) <...< Y n (t) < X n+1 (t), 

and the equations 

(8) Yi(t) = Vi + /3i(t A f ) + Lr(t At) — L+{t A f ), 

(9) Xi(t)=x i + ~ fi (tAt), 

where, 

f is the stopping time given by f — inf {t > : Xi(t) — Xi + \{t) for some i £ {1 , 2, . . . , n}}, 
/3i, /?2, • • ■ Pn, 7i, 72, • ■ ■ 7«+i are independent JVBrownian motions, 

the processes and L~ are continuous, non-decreasing and increase only when Yj = X i+ i and 
Y = Xi respectively, 



Lf(t) = [ 1{Y{ S ) = X i+1 (s))dLt(s) LT(t) = f l(Y(s) = Xi(a))dLi 
Jo Jo 



(s). 



Notice the difference between this process and [X, Y) is the reflection rule: here Y is pushed off X 
whereas it was X that was pushed off Y. 
Define a family (q";t > 0) via 



n+l,n 



(fO) q?((x,y),(x',y , ))=q?((x',y'),(x,y)) for (x,y), (x',y') e W 

The following proposition is proved by arguments exactly parallel to those just given in proof of 
Proposition [3 

Proposition 3. (q™;t > 0) are a family of transition densities for the process [X,Y\ killed at the 
instant f, that is to say, for t > and (x,y), (x',y') € W n+1 ' n , 

q^((x,y) 7 (x',y'))dx / dy' = Q^(X t G dx',Y t G dy';t < f). 

The duality, represented by l|l(J|) . between the transition semigroups of (X, Y) and (X, Y) is not 
unexpected, ft is consistent with general results, see for example DeBlassie 5 , and Harrison and 
Williams m, which show that, in a variety of contexts, the dual of a reflected Brownian motion 
is another reflected Brownian motion where the direction of reflection at the boundary is obtained 
by reflecting the original direction of reflection across the normal vector. This is precisely the 
relationship holding between (X, Y\ and here. 



3 An intertwining involving Dyson's Brownian motions 

It is known that Dyson's non-colliding Brownian motions can be obtained by means of a Doob h- 
transform. Let W n = {y G R™ : y\ < yi < . . . < y n }. Suppose that (Y t ;t > 0) is when governed by 
the probability measure P™ a standard Brownian motion in R n , relative to a filtration {Tt]t > 0}, 
starting from a point y G W n and stopped at the instant r = inf{< > : Yi(t) = Yj(t) for some i ^ 
j}. The transition probabilities of Y killed at the time r are given explicitly by the Karlin-McGregor 
formula, |12j . 



(11) 



P;(Y t edy';t<r) = P ?{y,y')dy' , 



for y,y' £ W n , where, with <p t again denoting the Gaussian kernel with variance t, 

(12) P?(y,y') = det{MVj - in); i < m < «}■ 

If the initial co-ordinates y satisfy y\ < j/2 < • • • < Vni then we may define a new probability measure 
by the absolute continuity relation 



(is) p r = i^T- p *" 

for t > 0, where /i„ is the function given by 

(14) h n (y)=H(y j - yi ). 



Under PJ/ + the process Y evolves as a Dyson non-colliding Brownian motion, that is to say r is 
almost surely infinite and the stochastic differential equations Jy hold. The transition probabilities 

(15) P n y ' + (Y t e dj/';t < r) = p? + (y, y')dy', 
are related to those for the killed process by an /i-transform 

(16) P?' + (y,y') = ^-Pny,y% 

h n {y) 

for y, y' € W n \ dW n . Finally we recall, see O'Connell and Yor, ^H], that we may describe Pq' + , the 
measure under which the non-colliding Brownian motion issues from the origin by specifying that it 
is Markovian with transition densities t > 0) and with the entrance law 

(17) PZ + (Y t edy)=tf(y)dy, 

for t > 0, given by 



?(») = ^~" 2/2 ex P | - E % 2 /(2*) } I Ufa ~ Vi) 



(18) £ 

with the normalizing constant being Z n = (27r)™/ 2 IIj<n i- 

Now suppose that [X , Yj is governed by the probability measure Q™ y defined in the previous 
section. Recall that > 0) are the transition densities of the process killed at the time r = 

inf{£ > : Yi(t) = Yj(t) for some i ^ j}. Suppose the initial co-ordinates y of Y satisfy y\ < yi < 
... < y n , then we may define a new probability measure Q"^ by the absolute continuity relation 

(19) q: . + = y^). Q: ,„ 

for t > 0. It follows from the fact that under Q™ y the process Y evolves as a Brownian motion stopped 
at the instant r, that h n (Y tAT ) is a martingale, and that this definition is hence consistent as t varies. 
Under the measure Q™'^~, the process Y now evolves as a non-colliding Brownian motion satisfying the 
stochastic differential equation JTJ, whilst the process X satisfies (QJ. The corresponding transition 
densities (<?™' + ; t > 0) are obtained from those for the killed process by the ^-transform 

(20) q? + {{x,y),{x',yi)) = } ^ q ?{{x,y),{x>,y>)) 
for (x, y), (a/, y') £ w n+1,n with the components of y all distinct. 

Lemma 4. The family of probability measures with densities given by {y™\t > 0) on W n+1 ' n , given 
by 



form an entrance law for (q™ , t > 0), that is to say, for s,t > 



<{x,y) = 



(w')= I v?(w)q?' + (w,w')dw. 



Accordingly we may define a probability measure Qo'g - , under which the process (X, Y\ is Marko- 
vian with transition densities > 0) and with the entrance law 

(21) Qo; + {X t G dx, Y t G dy) = v? (x, y)dxdy. 

It is easy to see that under this measure [X, Yj satisfies the equations (QJ and (QJ , starting from the 
origin x = 0, y = 0. Presumably any solution to and starting from the origin has the same 
law, but we do not prove this. 

We may now state the main result of this section. 

Proposition 5. Suppose the process (X t , Yt',t > 0) is governed by Qo'o" then the process (X t ; t > 0) 
is distributed as under Pq +1 ' + , that is as a Dyson non-colliding Brownian motion in W n+1 starting 
from the origin. 

This result is proved by means of a criterion described by Rogers and Pitman |18j for a function 
of a Markov process to be Markovian, see Carmona, Petit and Yor, for father examples of 
intertwinings. For x G W n+1 let W n (x) R" : Xi < y\ < . . . < y n < x n+ i, and define 

(22) A"(x,y) = n!- My) 



' h n+1 (x) ' 



for x G W n+1 \ dW n+1 and y G W n (x). The normalizing constant being chosen so that \ n (x, •) is 
the density of a probability measure on W n (x). This follows from the equality 

(23) / h n (y)dy = —h n+1 (x), 



which is easily verified by writing h n (y) = detjy^ x ;l < i,j < n}. The proof of Proposition [5] 
depends on the following intertwining relation between (q™' + ;t > 0) and > 0), for all 

t> 0, x G T^ n+1 \ dW n+1 , and (a;', y') G 

(24) / \ n (xM< + ((x,y) ) (x',y'))dy = p? +1 ' + (x,x')\ n (x',y'). 

JW n {x) 

This by be verified directly using the explicit formula for g™ given in the previous section. Alter- 
natively the following derivation is enlightening. Recall that if (X t ,Y t \t > 0) is governed by Q™ 

then the process (X t ;t > 0) is a Brownian motion stopped at the instant f = infji > 0;X; = 

Xj for some i ^ j } . Consequently the transition probabilities of the killed process satisfy 

(25) f q?({x,y),{x',y'))dy' =p^ +1 {x,x'). 

JW"(x') 

Now using the duality between g™ and g™ and the symmetry of p™ +1 we may re-write this as 

(26) f q?((x',y'),(x,y))dy' =p?+\x',x). 

JW n (x') 

Finally to obtain l)24|) we swop the roles of (x,y) and (x',y') and use the expressions for g™' + and 
^n+i,+ ag ^.^ rans f orins As a fi rg t application of the intertwining we have the following. 

Proof of Lemma^ Notice that v™{x,y) — fi™ +1 (x)\ n (x,y). Hence, by virtue of the intertwining 
and the fact that t > 0) is an entrance law for t > 0) we have, 

/ dxdy Vs(x 1 y)q?' + ((x,y), {x',y')) = 

J W n + l, n 

dxp n + l {x)( dyX n (x,y)qf' + ((x,y),(x',y')) = 

W n + 1 JW n (x) 

dx tf+\x)p? +1 ' + (x,x')\ n (x',y') = tf+\x')\ n (x',y') = v? +s (x',v')- 

□ 



A similar argument, following proves the proposition. 

Proof of Proposition^ For a sequence of times < t\ < t 2 < . . . < t n , repeated use of the inter- 
twining relation gives, 

Qo;a(x tl £A 1 ,X t2 £A 2 ,...,X tn £ A n ) = 

dxi... dx n dyi... dy n v\\ (x u yi)q^- tl Vi), ( x 2, 2/2)) • • • 

Ax J A n JW"(xx) JW n {x n ) 

■ ■ ■ C-t»-i (( x n-l($/n-l)> ( x n, Vn)) = 



dxi... dx n dyi... dy n pt" 1 +1 (xi)A"(xi, yi)q%± tl ((xi, yi), (x 2 , 2/2)) • ■ ■ 

Ai JA„ JW n (a;i) JW n (x n ) 

■ ■ ■ C-t»-i ((xn-i,y n -i), (x n , y n )) = 
I dx\...l dx n dy 2 ... dy n ^ 1 (x\)p^_^ {xi,x 2 )X l {x 2 ,y 2 ) . . . 

J A x J A n JW"(x 2 ) JW n {x n ) 

■ ■ ■ C-t»-i (( X n-l,y n -l), (Xn, Vn)) = 

dxi... \ dx n / dy n /i™ i +1 (a;i)p™ 2 t 1 t ' 1 + (xi, x 2 ) . . -Pt^-t^-x ( x n-i,x n )X n (x n , y n ) 

\ J A n JW n (x n ) 

dx x ... \ dx n ^™ +1 (xi)p™ 2 ^ 1 4 1 + (a;i, x 2 ) ■ ■ -Pt^-il-x ( x n-l,%n)- 
Ax JA n 

□ 

Notice that in the above proof, if we integrate y n over some smaller set than W n (x n ) we find 
that 

(27) Qoo + (^n G A\X tl ,X t2 ,...,X tn ) = f X n (X tn ,y)dy. 

J AnW"(X tn ) 

This may be interpreted as the following filtering property: the conditional distribution of Y t given 
(X; s < t) is given by the density X n (X t , •) on W n (X t ). 

4 Brownian motion in the Gelfand-Tsetlin cone 

Proposition [S] lends itself to an iterative procedure. Let K be the cone of points x = (x 1 , x 2 , . . . x n j 
with x k = (xf , x 2 , . . . , Xf".) £ R fc satisfying the inequalities 

(28) xt l <xf<x%l 

K is sometimes called the Gelfand-Tsetlin cone, and arises in representation theory. We will consider 
a process X(t) = (X 1 (t) , X 2 (t) , . . . X n (t)) taking values in K so that 

(29) X*(t) =x\ + iht) + £?'-(*) - (*), 

where (7^(^)5 1 > 0) for 1 < k < n, 1 < i < k are independent Brownian motions, and (L^' + (t); t > 0) 
and (^L^' 1 {t);t > 0) are continuous, increasing processes growing only when X^(t) = X^~ x {t) and 
X*(t) = X^Z\{t) respectively, the exceptional cases L^' + (t) and L x '~{£) being identically zero for all 
k. For initial co-ordinates satisfying x\ < x k +1 for all k and i, we may give a pathwise construction, 
as in Section 2, based on alternately using the Skorokhod construction to reflect Xf downwards from 
X\~ x and upwards from X)[Z\- The potential difficulty that X\~ x meets X*Z\ does not arise. 

In order to construct X starting from the origin we use a different method. First we note that if 
the pair of processes (X,Y), governed by the measure Qg'o^ satisfies equations {Q and 10}, then, 

(30) the Brownian motion 7 = (71, 72, ■ ■ ■ 7n+i) is independent of Y . 

By repeated application of Proposition |S1 there exists a process (X(t);t > 0), starting from the 
origin, such that 



the process (X k (t);t > 0) is distributed as Pq'~'~, for k — 1, 2, . . . , n, 

the pair of processes (X k+1 (t), X k (t); t > 0) are distributed as under Q 'o~j for k = 1, . . . ,n — 1, 

for fc = 2, . . . , n— 1 the process (X fe+1 (i); f > 0) is conditionally independent of (X 1 ^), . . . , X* -1 (t); 1 > 
0) given (X k {t);t > 0). 

By its very construction the process X satisfies the equations (|29|) . for some Brownian motions j k , 
which by the observation i|3U|) are independent. Even starting from the origin, pathwise uniqueness, 
and hence uniqueness in law hold for X. Consequently we may state the following proposition. 

Proposition 6. The process (X(t);i > 0), satisfying l|29|l . if started from the origin, satisfies for 
each k = 1, 2, . . . , n, 

(X {k \t)-t > 0) is distributed as under P^. 

The conditional distribution of (X k (t);t > 0) given (X k - 1 (t);t > 0) facto rizes in a Markovian 
fashion into the product of the conditional distribution of (X fe (s); < s < t) given (X k ~ 1 (s); < s < 
i), and the conditional distribution of [X k {u);t < uj given (X k - 1 (u);t < u) and X k (t). From this 
factorization we deduce that, for any t > 0, and k = 2, . . . , n — 1 the process (X fe+1 (s); < s < t) 
is conditionally independent of (X 1 (s), . . . X k ~ 1 (s); < s < t) given (X fe (s);0 < s < t). For 
any x k £ W k we will denote by K(a; fc ) the set of all (x 1 , x 2 , . . . x k ~ 1 ) such that for all i and j, 
%i + < ^ < • The — l)/2-dimensional volume of K(x fe ) is given by 

' h k (x k ). 



Ylj<k ■?• 



Recall from (I27|l that the conditional distribution of X k (t) given (X k+1 (s);0 < s < t) has the 
density X k (X k+1 (t) , •) on W k (X k+1 (t)). Combining this with the conditional independence property 
we deduce that the conditional distribution of (X *(<;), X 2 (t) . . . X k (t)) given (X k+1 {s)- < s < i) is 
uniform on K(JT fe+1 (t)). Finally using the fact that the distribution of X n (t) is given by the density 
/z" on W n we deduce that the distribution of X(t) has the density 



(31) m»(x) = (27r)-"/ 2 t- 2 / 2 exp - J](^) 2 /(2t) i J]^" " x ^ 




with respect to Lebesgue measure on K. Baryshnikov, 2 , studies this distribution in some detail. 
Let (H(t);t > 0) be a Brownian motion in the space of n x n Hermitian matrices, and consider the 
process (H 1 ^), H 2 (<),... H n (t); t > 0) where H k (t) is the k x k minor of H(t) = H n (t) obtain by 
deleting the last n — k rows and columns. It is a classical result that the eigenvalues of H k ~ 1 (t) are 
interlaced with those of H k (t). Baryshnikov shows that, at any fixed instant t > 0, the distribution 
of the eigenvalues of H l (t), H 2 (t), . . . H n {t) is given by the density (|31|l . However it is not the case 
that the eigenvalue process is distributed as the process (X(t);£ > 0). 

O'Connell, ^Bj) describes another process (T(t);t > 0) taking values in K which is constructed 
via certain explicit path transformations. This process arises as the scaling limit of the RSK corre- 
spondence. The process X described above has several features in common with T. For each fc, the 
subprocess (T k (t)- 1 t > 0) evolves as Dyson fc-tuple starting from zero. Additionally 

(32) (x 1 1 w,x 2 2 (i),...,x:(t);t>o)' i = t (ri(t),r 2 (i),...,r^i);t>o), 

but remarkably all other components Tf with I < k are given by explicit deterministic transformations 
applied to the processes T\, T\, . . . T™. A feature that X certainly does not share. 
Notice that, for k > 2, 

(33) X k (t)= 1 k (t) + L k k -(t), 

where L^ _ (t) grows only when X k (t) = X k z\{t). On applying the Skorokhod lemma, see Chapter 
VI of we find that 

(34) L^(t) = sup(A fe fe I 1 1 ( S )- 7 , fc ( S )). 



Iterating this relation we obtain 

k 

(35) X k k (t)= sup £{7i(ti)-7ft*<-i)}, 

Q=t <t 1 <t 2 <...<t k =t i=1 

which in the light of Proposition[fj] proves the identity ©. This is essentially the same argument for 
J2J as given by O'Connell and Yor, with Proposition replacing the corresponding statement 
about r. 

We close this section by noticing that (X{ (t) , Xf(t) , . . . ,X"(t);t > 0) is Markovian and giving 

(n) 

an explicit formula for it transition probabilities. For n > 1 let $J denote the nth order iterated 
integral of the Gaussian density fit , 

(36) i> t (n) (y)= f ^"^V ^dx, 



(n-1)! 



and for n > let $[ denote the nth order derivative of (fit- Define for x, x' <E W n , 



(37) r t {x,x') = det{$f~ J, (x; - Xi ); 1 < i, j < n}. 

Lemma 7. For any f : W n — > R which is bounded and continuous and zero in a neighbourhood of 
the boundary of W n , 



lim / r t (x,x')f(x')dx' = f(x), 

uniformly for all x S W n . 

Proposition 8. The process (Xf (t), Xf (t), • • • ,X™(t);i > 0) satisfying l|33|) is Markovian with 
transition densities given by r t (x,x'). 

Proof. For a fixed a:' € R, and any n, the function (i, x) i— > (x' ~ x) solves the heat equation on 
(0,oo) x R. From this we easily see that for a fixed x' £ R", that the function (t,x) i— > r t (x' — x) 
solves the heat equation on (0, oo) x R". Moreover if Xi = x.- L -\ for any i = 2, 3, . . . , n then the ith 
and [i — l)th rows of the determinant defining -^r-r t (x, x') are equal and hence this quantity is zero. 

Let / : W n — > R be a bounded, continuous and are in a neighbourhood of the boundary of W n . 
Then define a smooth function F on (0, oo) x W n via 



F(t,x) = [ r t {x,x')f{x')dx'. 



' w 

By virtue of the above observations regarding rt, and differentiating through the integral, we find 
that 

n B 2 F BF 
IE^M=a7/M on(0,oo)x^", 

i=l 1 

with the boundary conditions 
OF 

— — {t, x) = whenever Xi — Xi-i for some i = 2, 3, . . . , n. 

UXi 

Let X denote a process governed by a probability R^, with components Xi(t) < X<z(t) < . . . < X n {t) 
satisfying the equations Xuif) = x^ + jk(t) + L k (t), where 7^ are independent Brownian motions 
and L k is an increasing process growing only when Xk(t) = Xk-i(t), with L 1 being identically zero. 
Fix T, e > 0. Applying Ito's formula, we find that the process (Fir + e — t, X t ) ; t G [0, T]J is a local 
martingale, which being bounded is a true martingale. Thus 

F(T + e,(x))=R x [F(e,X(T))]. 

Appealing to the previous lemma, we may let e J, and so obtain, 

F(T,x)=R x [f(X(T))], 

which, since it is clear the distribution of X(T) does not charge the boundary of W n , proves the 
proposition. □ 



In view of Proposition we obtain from r t by a simple integration the following expression for 
the distribution function of the largest eigenvalue of H(t): 

(38) P^ + (X n <x)= det{$^' +1) (x)}. 
Possibly this can be checked directly using the Heine identity. 

5 Coalescing Brownian motions 

In this section we consider the joint distribution of a family of coalescing Brownian motions. Fix 
zi < z 2 < ■ ■ ■ < z n and consider the process of n coalescing Brownian Motions, 

t Z t = (Z t (zi), ■ ■ . Z t (z n )), 

where each process (Zt(zi); t > 0) is a Brownian motion (relative to some common filtration) starting 
from Zo(zi) — z%, with for each distinct pair i ^ j the process 

1 1-> ~^\ Zt ( z *) ~ Z t( z j)\ 

being a standard Brownian motion on the half-line [0, oo) with an absorbing barrier at 0. Thus 
informally [Zt{zi);t > 0) and (Zt(zj)',t > 0) evolve independently until they first meet, after which 
they coalesce and move together. Such families of coalescing Brownian motions have been well- 
studied, for some recent works concerning them see |5] and 

For a fixed t > 0, the distribution of Z t (z) is supported on W n . That part of the distribution 
supported on the boundary of W n corresponds to the event that coalescence has occurred. Whereas 
the restriction of the distribution to the interior W n (corresponding to no coalescence) is given by 
Karlin-McGregor formula : 

(39) P(X t (z t ) G dz[ for all i) = det{(^(^ - Zi)}dz'. 

In fact we can bootstrap from this result to a complete determination of the law of Z t (z), which can 
be expressed in the following neat way. 

Proposition 9. For z, z' £ W n , the probability 

P(Z t (z t ) < z[ forl<i<n) 

is given by the determinant of an n x n matrix with (i,j)th element given by 

Qtiz'j-Zi) ifi>j, 
$t(z'j-Zi)-l ifi<j, 

where 

* t (z)= f ^=exp{-y 2 /(2t)}. 

J -oo V ^t 

Proof. First we note that by integrating the Karlin-McGregor formula we obtain 

(40) P{Z t ( Zl ) <z[< X t (z 2 ) < z ' 2 <...< z' n _ x < Z t (z n ) < z' n ) - det{$ t (z; - z t )}. 

We are going to obtain the desired result by showing how the indicator function of the event of 
interest 

{Z t (zi) < z[, Z t (z 2 ) < 4> • • • j Z t (z n ) < z' n ) 
can be expanded in terms of the indicator functions of the events of the form 

{Z t {z, L{1) ) < z' j{x) < Z t {z l{2) ) < z' m < ... < 4( s -i) < Z t {z l(s) ) < z' j{s) }, 

for increasing subsequences of indices i(l), *(2), . . . , i(s) and j(l), j(2), . . . ,j(s). To this end I claim 
firstly that, whenever z, z' £ W n , 

(41) det{l(zj < z'j)} = l(zi <z' 1 <z 2 <z' 2 <...<z n < z' n ). 



I claim secondly that 



(42) det {^ z r*' j) Zi) \ ^} = < *i, *2 < 4, •••,*» < z' n ). 

To prove the first claim take the matrix M — {l(zi < z'j)}, and subtract from each column 
(other than the first) the values of the preceding column. The diagonal elements of this new matrix 
are 

l(zi < z[) - l(zi < zU) = 1(4-1 < ^ < z[)- 

adopting the convention that z' = — oo. Thus the product of these diagonal elements gives the desired 
result. We have to check that in the expansion of the determinant this is the only contribution. 
Suppose that p is a permutation, not the identity. Then we can find i < j with p(i) > i and 
p(j) < t- Consider the product of the (i,p(i))th and (j,p(j))th elements of the matrix (after the 
column operations). We obtain 

H Z p[i)-l < Z i< Zp^Hz'pU)-! < Z 3 ^ Z ' P (3))- 

This can only be non-zero if both 4(i)-i ^ Zi anc ^ z i — z 'p(jY kut Zi — z i so ^ s wou ld imply 

z 'p(i)-i z 'p(j)' ^ n view °f the f ac t P(fy ~ 1 — PC?) this ^ s impossible. 

Consider the matrix N appearing in the second claim. The product of its diagonal elements 
gives the desired result. To show that this is the only contribution to the determinant, take p a 
permutation, not equal to the identity and i < j with p(i) > i and p(j) < i, as before. Then the 
product of the (i, p(i))th and (j, p(j))th elements of the matrix is 

_1 < ^(zj < z' p{j) ) 

Since Zi < Zj for this to be non-zero we would have to have z', ^ < z' p ,^ , which is impossible for 
P{i) > PU)- 

Let T = { — l(j > i)} be the upper triangular matrix so that N = M+T and consider the Laplace 
expansion of det(M + T) in terms of minors. For increasing vectors of subscripts i and j let M[i, j] 
denote the corresponding minor of M and let T[i, j] be the complementary minor of T so that 

det(JV) = det(M + T) = ^(-l) s(iJ) M[i,j]T[i,j], 

ij 

for appropriate signs s(i,j). Evaluating det(iV) via the second claim, and the minors M[i,j] via 
(general versions of ) the first claim we have obtained an expansion of l(zi < z[ , Z2 < z' 2 , . . . , z n < z' n ) 
as a linear combination of of terms of the form l(zj(i) < ^Jm < z «( 2 ) — 4(2) < ■ ■ • — 4(«4 

To complete the proof replace, in the above expansion, Zi by Z t {z{) and take expectations. On 
the lefthandside we obtain P(Z t (zi) < z[, Z t (z 2 ) < z' 2 , ■ ■ ■ , Z t (z n ) < z' n ). On the righthandside we 
have a linear combination of probabilities: P{Zt(zi(x)) < z'^^ < Z t {zi( 2 ^) < z'^ 2 s < . . . < 4(s-i) < 
Zt(zi( s )) < z 'j( s )) eacri °f which can re- written by means of the integrated Karlin-McGregor formula 
as a minor of the determinant det{$(z^ — z;)}. And to finish we notice that the righthandside is now 
the Laplace expansion of the determinant of the sum of matrices {^{z'j — z^} and { — l(j > i)}. □ 

The expression just obtained for the distribution of coalescing Brownian motions is closely related 
to the formula for the transition density of the interlaced Brownian motions given by Proposition 
In fact it is easily verified that 

an »+l 

(43) q?{(x,y),(x',y>)) = (-1)" — - — ^ - ■ ■ ^- P(Z t ( Xi ) < x^Z t { Vj ) < y) for alH, j). 

This represents a duality between the the interlaced Brownian motions and coalescing Brownian 
motions which generalizes the well-known duality between Brownian motion on the half- line [0, oo) 
with a reflecting Barrier at zero, and Brownian motion on the half-line with an absorbing barrier at 
zero. 

There is interesting alternative way of expressing the equality 143(1 . The Arratia flow or Brownian 
web is a infinite family of coalescing Brownian motions, with a path starting from every point in 



space-time. Let t € [s, oo) i— > Z s j(x) denote the path starting from (s, x). It is possible to define on 
the same probability space a dual flow with paths running backwards in time: s S (—00, t] 1— > Z s j(x) 
being the path beginning at (t,x). For the details of this construction see [21] and [HJ. The flow Z 
and its dual Z are such that for any s,t,x and y, the two events Z s j(x) < y and Z St t(y) > x differ 
by a set of zero probability. Using this we may rewrite 143J1 as 

(44) <ff({x,y), (x',y'))dx'dy = P(Z Q , t (xi) e dx'^Z^y)) G d% for all i, j). 

This seems to fit with the fact that the paths of Z are "reflected off" those of Z, see ^5] an d EH|- 

6 Proofs of two lemmas 

Proof of Lemma^ The contribution to the determinant defining r t (x, x') coming from the principal 
diagonal is equal to the standard heat kernel in R™. The lemma will follow if we can show all other 
contributions to the determinant are uniformly negligible as t tends down to 0. Choose e > so that 
the function / is zero in an 2e-neighbourhood of the boundary of W n . Then consider a contribution 
to the determinant corresponding to some permutation p which is not the identity. There exist i < j 
with p(i) > i and p(j) < i, and the contribution corresponding to p consequently contains factors of 
<$>¥~ p{t)) (x' p{i) - x t ) and (2/ - Xj). Noting that j - p(j) > and i - p(i) < we see that 

on the set {x' p u\ — Xi > e} U {x' p /j\ ~ Xj < — e} at least one of these factors, and indeed the entire 
contribution, tends to zero uniformly as t tends down to zero. But on the complement of this set we 
have x',^ <Xi + e<Xj + e< x' p ,js + 2e, and p(j) < p(i) implies that x' p ,^ < x',^, so we see that 
x' is within the 2e-neighbourhood of the boundary of W n , and does not belong to the support of /. 
This proves the lemma. □ 

Proof of Lemma^ It is convenient to write z\ — Xi,Z3 = X2, ■ ■ ■ Z2n+i = x n , and z-i — 2/1,24 = 
J/2, ■ ■ ■ j Z2n — Vnt with a corresponding change of notation for x\ and y[ also. Now reorder the columns 
and rows of the determinant defining q" so that the (i, j)th entry is a function of the difference z'j — Zi. 
We may now argue in the same way as in the preceding proof. Choose e > so that the function / is 
zero in an 2e-neighbourhood of the boundary of W n+1 ' n . Consider a contribution to the determinant 
corresponding to some permutation p which is not the identity. There exist i < j with p(i) > i and 
p(j) < h and the contribution corresponding to p consequently contains factors which are functions 
of z'p(f) — Zi and z' p ,-\ — Zj. Noting that j — p(j) > and i — p(i) < 0, and checking the entries of the 
determinant above and below the diagonal we see that on the set {2^ — z% > e} U { z ' p (j) — z j < — e } 
at least one of these factors, and indeed the entire contribution, tends to zero uniformly as t tends 
down to zero. As above, this proves the lemma. □ 
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